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We calculate the harmonic-gauge even / = 1 mode of the linear metric 

perturbation (MP) produced by a particle in a weak-field circular orbit around 

a Schwarzschild black hole (BH). We focus on the Newtonian limit, i.e. the 

limit in which the mass M of the central BH approaches zero (while fixing 

the orbital radius and the small-object mass), and obtain explicit expressions 

for the MP in this limit. We find that the MP are anomalous in this limit, 

namely, they do not approach their standard. Coulomb-like, flat-space values. 

Instead, the MP grows on approaching the BH, and this growth becomes 

worse as M decreases. This anomalous behavior leads to some pathologies 

which we briefly discuss. We also derive here the next-order correction (in 

the orbital frequency Q) to the MP. 



I. INTRODUCTION 

One of the most useful gauges for analyzing linear metric perturbations (MP) h^^ over 
a given background metric g^v is the harmonic gauge. This gauge may be regarded as the 
gravitational analog of the electromagnetic Lorentz gauge. The linearized Einstein equa- 
tions in this gauge reduce to a hyperbolic system [Eq. (2) below], which allows an elegant 
formulation of the initial-value problem. Also, the phenomenon of gravitational self force 
[1] [2] was originally formulated within this gauge. In this paper we shall show that the 
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harmonic gauge behaves in an unusual manner when a source is present in the weak-field 
region of a Schwarzschild black hole (BH). [3] 

The harmonic gauge is defined by the condition 

= 0, (1) 

where h^j, is the trace-reversed MP, 

The linearized Einstein equations then reduce to 

h^,.;^ + 2R^^p,h'^^ = -167r7;., (2) 

where Ranpv is the Riemann tensor associated with the background metric g^i,. 

We shall consider the gravitational perturbation produced by a point particle of mass 
11. The simplest such problem is that of a static particle in flat space. In this case the 
trace-reversed harmonic-gauge MP takes the simple form 

^« = T' 

where d is the distance from the particle (in the latter's rest frame), with all other compo- 
nents of h^i, vanishing. This yields 

hij.li = ~^Sfj^i,, (3) 

(in Cartesian coordinates), where 5^,^ is 1 for jj, = u and zero otherwise. 

Suppose now that the background geometry includes a BH of mass M, located at a 
distance R from the particle. Obviously the fiat-space solution (3) will no longer be vahd. 
In this paper we shall address the following question: What will happen to the harmonic- 
gauge functions h^i,{r") at the Newtonian limit where M is taken to be arbitrarily small 
while R and the evaluation point are kept fixed? Will they approach the fiat-space 
solution (3)? We shall analyze this problem (specifically for circular orbits in Schwarzschild 
spacetime), and show they do not. 



The same question may be asked for various kinds of perturbation fields in a weak-field 
BH spacetime, and also to various modes of such fields. Assume, for example, that the 
point particle is electrically charged. Will the electromagnetic field produced by the particle 
approach its simple Coulombic flat-space form, in the presence of a distant BH, in the 
Newtonian limit M — > 0? We shall say that a field "behaves normally" in the Newtonian 
limit if it does approach its flat-space value as M — > 0. If, on the other hand, the field 
(or a particular mode of the field) approaches a different limit as M — > 0, we shall refer 
to this field (or mode) as anomalous. In this paper we shall show that the polar I — 1 
mode of the harmonic-gauge MP is anomalous; and so is the overall harmonic-gauge MP 
(i.e. the sum over all modes). As a consequence, when a small BH is located at a large 
distance R from a finite-mass source, the harmonic-gauge MP is pathological in some sense, 
and should be used with great care. In particular, if /iR » (even though fi « M ) 
the harmonic-gauge MP grows large (i.e. hfj,,, » 1) in the neighborhood of the BH, and 
the linear perturbation analysis actually breaks down. (Obviously, it is possible to use the 
linear perturbation analysis for the MP, even in this case — but not in the harmonic gauge!) 

On physical grounds, one may expect that all locally-measurable quantities will behave 
in a normal manner as M — > 0: It is unlikely that a local measurement will be significantly 
affected by the presence of a distant BH if the latter's mass is arbitrarily small. Therefore, 
the anomalous behavior of the harmonic-gauge MP should be attributed to a bad choice 
of gauge. That is, the harmonic gauge is advantageous in many cases, but it becomes 
problematic (or even invalid, as in the case fiR » M^) in some weak- field situations. 

The motivation for this work came from a recent analysis by Pfenning and Poisson [4]. 
They calculated the gravitational self force (as well as the electromagnetic and scalar self- 
forces) acting on a point particle in a weak-field curved spacetime. Their analysis, which was 
based on the harmonic gauge, indicated the crucial role of the I — 1 mode in understanding 
the (weak-field) gravitational self force, and hinted that this mode might behave anomalously 
in the weak-field limit. Here we confirm that this is indeed the case. 

Let us now formulate our problem in a more precise form. We shall consider a small 
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particle of mass // moving on a circular geodesic orbit around a Schwarzschild BH of mass 
M » ji. We use the Schwarzschild coordinates for the background, 

ds"^ = -(1 - 2M/r)de + (1 - 2M/r)-^dP + r'^{de'^ + sin^ Qdip'^). 

The circular orbit is assumed to be equatorial, i.e. 9 — 7r/2, and its radius is r = i?. The 
two constants of motion are denoted E = —ut and L = u^. The orbital frequency is 

n = ^^ {M/R^Y'\ 

CLv 

We shall consider a weak- field situation, i.e. R» M, and analyze the even I — 1 mode of 
the harmonic-gauge h^j,^. Particularly we shall calculate the limit M ^ (for fixed radius 
R and location r") of this mode. We shall show that this limit does not agree with the 
corresponding flat-space expression; that is, the even I — 1 mode of the harmonic-gauge MP 
is anomalous. 

We then go one step forward and calculate the first-order correction to this mode in the 
small parameter i?Q = {M/Rf/'^. 

A. The general method 

In our circular-orbit problem the energy-momentum source is restricted to a delta- 
function distribution at r = i?. At both sides r < R and r > R the MP is source-less. 
We denote the harmonic-gauge vacuum MP at r < i? and r > Rhy and respec- 
tively. 

To construct the desired harmonic-gauge functions h^^^ we proceed as follows. Accord- 
ing to an analysis by Wald [5], concerning the Z = 0, 1 modes of the gravitational pertur- 
bation over the Schwarzschild geometry, if in some range Ri < r < R2 oi the perturbed 
Schwarzschild background the gravitational perturbation is sourceless, then a gauge exists 
in which the even I — 1 part of the MP vanishes throughout this range. (In fact Wald's 
analysis is more general, and it addresses the Kerr case as well.) Thus, inside this range the 
even I — 1 MP may be described in terms of some even I — 1 gauge vector through 
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(4) 



If we further impose the condition 



- ;A 



(5) 



then the harmonic gauge condition (1) is satisfied by /i^j,. The wave equation (2) will then be 
satisfied too, because the derived h^u is just a pure-gauge perturbation. Thus, the problem 
of constructing the even I — 1 vacuum harmonic-gauge MP reduces to that of finding the 
even I = 1 solutions to Eq. (5). We shall refer to such solutions of Eq. (5) as "homogeneous 
gauge vectors" . 

Applying this procedure to each of the ranges r < R and r > R, we conclude that the 
two vacuum harmonic MP solutions h^^^ and h^^^ may be derived through Eq. (4) from two 
even / = 1 gauge vectors that solve Eq. (5) (one in each range), which we denote ^''^^ and 
respectively. 

We emphasize that there does not exist a global gauge vector ^iaiohai) ^j^g entire domain 
r > 2M such that h^y is globally derived from it via Eq. (4). In other words, the even / = 1 
metric is trivial (i.e. may be annihilated by a gauge transformation) locally va. r < R and 
separately in r > i?, but not in the entire domain r > 2M. 

In principle the construction of the harmonic MP consists of two main stages: (I) Con- 
structing the general homogeneous gauge vectors, namely the general even I — 1 solutions 
to Eq. (5); and (II) Finding ^^"^ and ^^^\ i.e. the specific, actual gauge vectors r < R 
and r > R, respectively. The choice of the actual gauge vectors ^j^"* should be made, in 
principle, by first constructing the general vacuum solutions h^^^ from the general homo- 
geneous gauge- vector solution at both sides, and then imposing the following conditions on 
h^i^^: (a) Appropriate asymptotic behavior at the event horizon (EH) and at infinity, (b) 
continuity at r = i?, and (c) A proper discontinuity of the r -derivatives of h^^^ a,t r — R, as 
determined from the energy-momentum distribution there through the field equation (2). 

In practice, the condition (a) is most easily implemented at the level of ^. As can be 
easily verified, if a gauge field is regular at the EH, the MP constructed from it via Eq. 
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(4) is regular at the EH too. The situation at the hmit r — > oo is similar: If a gauge field 
vanishes at infinity, the corresponding MP will vanish at this limit too. 

To the best of our knowledge the exact explicit solutions to Eq. (5) (for M 7^ and 
frequencies cu 0) are not known. Here, however, we are only interested in the small-M 
limit, i.e. the situation in which M « R. We shall thus calculate the MP at the leading 
order (M/R)'^, (Rfl)^, and also at the next order {Rfiy. To this end we shall only need 
an approximate expression for the gauge vector, accurate to a certain order in the small 
parameters M/R and Such an approximate expression can easily be constructed. 

In section II we shall construct the general homogeneous even / = 1 gauge vector, at the 
leading order [i.e. order (M/i?)°, (i?Q)°]. Then in section HI we correct this homogeneous 
solution at orders {M/RY and {RilY (these turn out to be the orders required for our anal- 
ysis). In section IV we consider the energy-momentum source associated with the massive 
particle and construct its even I — 1 mode. Then in section V we use this source, combined 
with the general homogeneous solutions derived in section HI, to construct the actual gauge 
vector at the two vacuum regions r < R and r > R. From this gauge vector we then 
derive in section VI the harmonic-gauge even I — 1 MP functions h^^^ — first at the leading 
order {M/R)° — (i?Q)°, and then at the next order (RQy. Then in section VI we explore 
the overall harmonic-gauge MP, i.e. the sum over the contributions of all modes, at the 
hmit M — > 0. We thus express the overall MP as the sum of the fiat-space term h^^^ and 
an anomalous term A/i^j, (independent of M). Finally, in section VIII we briefiy discuss 
the results. In particular we show that the harmonic-gauge MP becomes >> 1 (and hence 
invalid as a linear perturbation) if the particle is located too far from the BH. In Appendix 
A we present five of the six exact even, I — 1, static homogeneous solutions to Eq. (5) in 
the Schwarzschild background (the sixth solution, which we denote is not presented 

as it is too long). 
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II. HOMOGENEOUS GAUGE- VECTOR SOLUTIONS: LEADING ORDER 

A. Flat-space static solutions 

We start by exploring the static even I — 1 solutions of Eq. (5) in flat space. For each 
azimuthal number m {m gets the three values 1 < m < 1) there are six such solutions, which 
are naturally divided into three pairs. We denote them where the "+" 

label refers to solutions that vanish at infinity, and the "-" label refers to solutions that are 
regular at the origin. 

We first present the six m = solutions, which are axially-symmetric. The first pair is 
the "scalar mode", 

where 

0(-) _ fQQs0 ^ 0(+) — r~^cos^. 

[Note that 0^^^ are just the two axially-symmetric I — l,m — solutions to the scalar wave 
equation.] The next pair is the "electric mode" 

^p-) = r cos ^ , CP^^ = cos 6, (6) 

The last pair is 

ef-^^r^cos^, =r-^cos^, 
^f-)^lr'sme, ef+^ = -sine, 

= = 0- 
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The related m — ±1 modes are obtained by a straightforward rotation of the m — 
solution. For our purpose it is convenient to consider the real modes (i.e. sin and cos 
functions of (p rather than complex exponents). In one such combination of the m = ±1 
azimuthal numbers (the combination which is mostly relevant to the analysis below) the six 
solutions are as follows. The "scalar mode" is 

with 

= r sin ^ cos ip , 0^"*"^ = sin 9 cos ip, (8) 

the "electric mode" is ^ 

^ — rsinOsimp , ^P^'* = sin 6 sin (p (9) 

and the third pair is 

^(^-) = sin Ocosip , = sin 9 cos (11) 

^ — —-r^ cos9 cosip , ^g^^'' = cos 6* cos 9?, (12) 

if"'^ sin ^ sin , +^ = - sin ^ sin (13) 

^) = 0. (14) 
The other set of six m = ±1 solutions is now obtained by replacing tp ^ tp + n /2. 



^Here and throughout this paper we take this specific phase (sin (p as opposed to the phase factor 
cos if used above for the scalar mode) for the "electric mode" , because it is this phase which emerges 
from the matching to the source at r = i2, (see section V). 



B. Slowly-rotating weak-field solutions 



We now return to the problem of weak-field circular orbits in Schwarzschild. Our goal 
at this stage is to construct the relevant homogeneous solutions, at the leading order in M 
and fl, namely at orders {M/ Rf and At this order, all radial functions simply take 

their above static flat-space form. Since the orbit is equatorial, the m = solutions turn 
out to be irrelevant (because the even I = 1 part of the source term does not include m = 
modes — see section IV). The m = ±1 solutions for the slowly- rotating modes are obtained 
from Eqs. (7-14) by simply replacing (/?—>(/? — Q,t. The set of six leading-order solutions is 
again divided into three pairs. The "scalar mode" is 

= ' (15) 

with 

= rsinecos((^ - Vtt) , 0^+^ = r~^sinecos(^ - Vtt). (16) 
The "electric mode" is 

= r sin ^ sin((^ - VLt) , ^^^^ = r"^ sin 9 sin((^ - VLt) (17) 

with 

The last pair is 

^f-) = sin 9 cos((/? - Vtt) , = r"^ sin 9 cos((^ - Ot), (19) 

= -]-r^ cos e cos((^ - ^It) , = cos 9 cos((^ - Oi), (20) 

^(p-) = ir^ sin 9 sm{^ - Vtt) , ^f^^ = - sin sin((^ - Qt), (21) 

= 0. (22) 



As before, the other six m — ±1 solutions are obtained by transforming (p — > (/? + 7r/2. 
(However, this second set of m = ±1 solutions will not be involved in the construction 
below, for the same reason as for the set of m = solutions.) 

Note that we do not expand the functions cos(</7 — ^t) or sin((/7 — Qi) in Q, for the 
following obvious reason: Although Q is assumed to be small (compared to e.g. 1/R and 
1/r ), fit may be arbitrarily large. Only the radial functions are expanded in fl — and in 
this section we consider the leading order, which is independent of Q (or M). 

Our starting point in constructing the above leading-order solutions was the set of fiat- 
space static solutions. Alternatively, one may start from the exact static solutions on the 
Schwarzschild background, given in Appendix A, and then take the limit M — > 0. This 
leads to the same results, Eqs. (15-22). In the static Schwarzschild case, too, there are 
three "-" solutions that are regular at the EH, and three " +" solutions that vanish at 
infinity. It is important to note that at the hmit M — > the Schwarzschild's and "-" 
solutions properly approach the corresponding fiat-space and "-" solutions. That is, 
the Schwarzschild's / = 1 gauge fields are non-anomalous at the limit M — > (this is to be 
contrasted with the corresponding MP functions, which are anomalous, as we show below). 

III. HOMOGENEOUS GAUGE- VECTOR SOLUTIONS: HIGHER-ORDER 

CORRECTIONS 

In the construction of the actual gauge- vector solution, some of the above homogeneous 
solutions must be calculated to higher order in the small parameters M/r and rQ. The basic 
reason is the following: As we show below [see e.g. Eq. (35)], in constructing the actual 
gauge- vector solution the homogeneous solution enters with a prefactor propor- 
tional to {M/R)~^; and similarly, enters with a prefactor proportional to (RQ)~^. [All 
other homogeneous solutions enter with a prefactor proportional to (Rfl)'^ or even smaller]. 
Note the following anomalous situation: Whereas the harmonic MP scales at the leading 
order as (M/it!)°, the actual gauge vector scales as {M/R)~^. In order to obtain the MP 
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at the leading-order (M/R)^, we must determine the actual gauge vector at the same 
order. This requires us to determine C'^'^ with accuracy {M/rY (rf2)^. 

Furthermore, we would like to go one step forward and to evaluate the correction to the 
MP proportional to (RflY- This requires one to determine with the same accuracy {RQY- 
In turn, this requires us to determine Cjf'^ order (rf2)^ [this also includes corrections of 
the type {rQ){M/r) = (Mil)], at order (rO)^ and (M/r), and all other homogeneous 

solutions at order (rQY- 

In practice, the corrections to the radial functions of the homogeneous gauge vector only 
emerge in even powers of rfl (but any power of M/r ), and no odd powers of rfl are involved 
(this law brakes at some order in RD, which is beyond our present goal). This implies the 
following: 

To calculate the MP at order {RilY, we need to evaluate (the generating function 
for ^''^~^) at order (rJl)^ and {M/rY, ^"<i other homogeneous solutions at the leading 
order (rQ)° and {M/r)^. And in order to get one extra order (i.e. order RQ) in the MP, all 
we need is to calculate , too, at order {rVLY and {M/rY- 

Summarizing these considerations, since our ultimate goal is to calculate the MP at 
accuracy {RQYi we need to calculate 0*^'^ and ^^^"^ at order (rQ)^ and {M/rY, ^^<i 
other inhomogeneous solutions at the leading order {rVlY and {M/rY ■ 

Note that up to the order required here [i.e. {rflY {M/rY] there is no couphng 
between the effects of rotation (i.e. 7^ 0) and curvature (i.e. M 7^ 0). These effects first 
couple at order (rf2)(M/r) = MQ = (ril)^, which is not required in our analysis. This 
situation allows the following simple procedure. The 0{rflY corrections are obtained by 
exploring the flat-space Q 7^ radial functions; And the 0{M/r) corrections are determined 
from the exact static homogeneous radial functions in the M 7^ Schwarzschild case. Note 
that both kinds of corrections are only required for the two homogeneous solutions ^^"^"^ 
(or, in practice, (t)^~^) and In what follows we shall describe in more detail these two 

types of corrections. 
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A. Rotation-induced corrections 



Consider first the scalar solution ^("■^ We need to explore here the homogeneous solutions 
to the scalar / — 1 radial equation, for a mode cu = Q [specifically we consider here the 
solution with the angular dependence given in Eq. (16) above]. We only need here the "-" 
solution, i.e. the one regular at the origin. This solution takes the simple exact form 

sin(rr2) — rflcos{rfl) 



{rny 



r sin 9 cos(</? — Qt) 



(the factor 3 was introduced here in order for 0^"^ to comply with the above leading-order 
expression for In our analysis we shall only need the corrections to the radial function 

up to order (rfl)^, which read 

(r(])2" 



0(- 
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rsm9cos(ip — Qt). 



Next consider the exact solution ^jf ^ for a flat space but with Q 7^ 0. The equation 
imposed on is 

□ef -) = 0. 

In fiat space, each Cartesian component satisfies this box equation separately. Therefore, 
an exact solution to this equation is 

"sin(/-r)) - /•r)siii(/-r))' 



{rfty 



rsin6'sin(<^ — Q.t) . 



(23) 



We adopt this solution as the exact version of the leading-order expressions (17-18) for 
The decomposition of the radial function in (rfl) then yields 



^In the static flat-space case coincides with the I = 1 mode of the Lorentz-gauge electromag- 

netic four-potential (hence the phrase "electric mode"). In the $7 7^ case, as defined here 

no longer agrees with the Lorentz-gauge electromagnetic four-potential, as it is not divergence- 
free. Nevertheless, Eq. (23) is a legitimate extension of Eqs. (9-10) to the 7^ case, which 
considerably simplifies the calculations. 
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1 - 
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rsm9sm{ip — fit). 



B. Curvature- induced corrections 

To obtain the M/r corrections we explore the exact solutions for the homogeneous, static, 
I — 1 gauge vector in Schwarzschild spacetime (see Appendix A). For the scalar mode 0^"^ 
we have the simple exact solution 

0(-) — (j, _ ]\^^ gjj^ Q COs{(p — (fo) . 

For ^jf'^ we have the exact solution 



^^«-) = 5* (r - 2M) sin^sin(v9 - ifo) ■ 



C. The total corrections 



Summing the above corrections to the radial functions of and which emerge 

from both the rotation and the nonvanishing of M (and expanding them in the small pa- 
rameters rQ and M/r), we obtain 



1-^ 

r 
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r sin 9 cos{ip — Qt) 



(24) 



and 



(B-) 



1 - 



2M {rny^ 
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r sin6'sin(9? — fit). 



(25) 



These corrected expressions for (f>^~^ and will be used in section V, together with the 
leading-order expressions (15-22) for 0^"*"-', and ^<^'"^\ in the construction of the actual 
gauge vector 
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IV. THE ENERGY-MOMENTUM SOURCE 



To construct the actual gauge vector at both sides of r — R we need the I — 1 part of 
the energy-momentum distribution at r = R. 

The particle's equatorial circular orbit is given by 

r = i?, e = 7r/2, if = -u^r, t = u^t. 

For this kind of orbit u'^ and are constants that depend on M and R, and satisfy fl — u'^/u^ 
(without loss of generality we have set i = (/? = at r = 0). The energy-momentum tensor 
is 

= H{R^ sin euY^5{r - R)5{e - 7r/2)5{ip - nt)u^u^. (26) 

This tensor is now expanded in tensorial harmonics Y^l)'''^. Throughout this paper we use 
the notation introduced in Ref. [6] for this decomposition. For a general tensor 7)ji,(x") the 
tensorial mode decomposition takes the form 

oo I 10 
(=0 m=—l i=l 

The "coefficients" T*^*)'™ are given by 

r«'-(r,t) = ^«|r;^V''rir(^'<^)]*^M'^(^^'^'^)«i^^^^^^' (27) 

where 

77^" = (iia^[-l,l,r"^(rsin^)-2], 

and /c*^*^ is either +1 or —1. For a general / the harmonics i = 1, ...,7 are even and i = 8, 10 
are odd, but for Z = 1 the harmonic i — 7 vanishes identically. Since Ur — 0, the source term 
(26) does not have any T^j, component, and consequently the harmonics i — 2,3,5 vanish 
(see Eq. (13) in Ref. [6]). We are thus left with three even harmonics, i = 1,4,6: 
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± 

rril=l,even _ 

m=— 1 

= ^(l)+^(4)+^(6)_ 

(HI/ ' /Hi/ ' /Hi/ 

(Note that A; = 1 for i = 1, 6 and A; = — 1 for i = 4.) We shall now proceed to calculate the 
contributions of these three harmonics. 
The harmonic i = 1 is 

YA(l)Zm _ VA^mrf rt 

and a straightforward calculation yields 

T«'"^ = l_iE\R\y^6{r - R)[Y^"'{e = 7r/2,ip = Qt)]* . 

Since Y^'^ oc cos^ vanishes at ^ = 7r/2, the only relevant spherical harmonics are m = ±1, 
given by ^ 

Y^'^\e,(p) = Ti^3/87rsin^e^^^. 

One finds 

ya)i,±iT(i)i'±i = (3/87r)//£;2(i?2^*)-i5(r _ i?) sin ^e±'(^-"*) 5^(5* , 

hence 

Z^ij = {3/47r)i^E\R'^u')-^S{r - R) sin^ cos((^ - . 
The i = 4 tensor harmonic is 

^2l[l + 1) 

In calculating T^'^^"* through Eq. (27), the term oc F^,"* of Yjf^^'^ multiphes Tte, which 
identically vanishes for an equatorial orbit. The only contribution to T*^')'"* thus comes from 



^We adopt here the notation used in Ref. [7] for the scalar spherical harmonics. 
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the term proportional to The latter vanishes for m — 0, hence again there will only be 



contributions from m — ±1. The explicit form of 1^^^^'"* for these modes is 

X [± cos eiS'X + ^X) + ' ^(^I'^T + '^T'^t)]- 



y(4)i,±i 



After a straightforward calculation one finds 
which yields 

cosesin(^ - n()(5'X + S'X) + sinflcos(^ - fitXi^iJ + i^i^) 
Finally, the i = 6 harmonic is 

Y^t""^ = ^Y^^iSX + sin^ eS^^S^). 

Again, T^^^^'^^ will not get a contribution from the term oc in rj6)/m^ because Tge 
vanishes. The contribution from the other term oc 5^5^ will be proportional to T^^, which 
in turn is proportional to oc Q^. Since our analysis is limited to the order Q^, it will not 
include contributions from i = 6. 

In summary, at the desired accuracy (fi^), the energy- momentum distribution is 

^^=l,e.en^^a) + ^(4). (28) 

Since we ignore here 0{M/R) correction terms, we may substitute li* = £■ = 1 and B?Q,, 
yielding 

Z« = (3/47r)/ii?-25(r - R) sin 9 cos{^ - Vtt)5l5l (29) 

and 

Zf} = -{3/87r)i^nS{r - R) x 

cos e sm{ip - Qi) (5* 5^ + 5l5l) + sin 9 cos{ip - Qt) (5^5^ + 5'^5l) 

Note that Zj^J is of the leading order 0° and Zj^J is O(O^). 

16 



V. THE ACTUAL GAUGE VECTOR 



In both ranges r < R and r > R the perturbation is sourceless, hence in each of these 
ranges the even / = 1 mode may be described by a certain homogeneous gauge vector [i.e. 
a certain even I — 1 solution to the equation (5)]. We denote this actual gauge vector by 
a,t r < R and a,t r < R. Each of these vector fields may be expressed as a certain 
combination of the basis homogeneous solutions constructed in sections II and III. In order 
for to vanish at r — > oo, it must be composed of the three "+" solutions solely. Similarly, 
in order for to be regular at the EH, it must be composed of the three "-" solutions 
solely. ^ Thus, ^^"^^ must take the form 

- a^^^ei^^^ + fc^^^er^ + ^^^^er^- (30) 
The MP functions h^^^ at the two sides oi r — R are constructed from ^j^^ via 

/^2^ = -(ei?+eS)- (31) 

The six coefficients a'^^\b^^\d^^ are to be obtained by imposing two matching conditions 
on the MP: 

(i) hf^u must be continuous at r = i?: 

(ii) The radial derivatives h^i^^r will undergo a certain jump at r = i?. This jump is 
conveniently described in terms of the trace-reversed MP, /i^,^, through the quantities 

^Recall that the three approximate "-" solutions constructed above may be obtained from 
the corresponding three static "-" solutions in Schwarzschild [by taking the leading order (s) in 
M/r, and adding the rotation-induced correction terms when required]; and the three "-" static 
Schwarzschild's solutions are regular at the EH by construction (see Appendix A). 
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The jump is to be determined from the even I — 1 energy-momentum distribution at r — R, 
Eq. (28). Applying the linearized Einstein equation (2) to this thin shell (approximating 
grr ^ one finds the required jump conditions: 



< htt,r >= -12^R~^ sin e cos{(p - ilt) (32) 



at order and 



< hte^r >—< hot,r >— 6 lift cos 9 sm{(p — fit), (33) 

< hf^^r >=< h^t,r >= 6liflsm9 cos{(p — fit) (34) 

at order fl^, with all other components of < hf^^^r > vanishing. 

The calculation of the six coefficients a^^-* , 6*^^^ , c*^^^ thus proceeds as follows. One first ap- 
plies Eq. (31) to the actual gauge vector written in the form (30), then apply the continuity 
condition (i), and then the jump condition (ii). The calculation is lengthy but straightfor- 
ward, and we shall skip here the details. The values of the above six coefficients are found 
to be 
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^^''^ = T7^f^^ (35) 



and 



a(+) = 0, = -^flfi , cW = -2i?/i . (36) 
5 

The actual gauge vector at r < it! then becomes 

^* ^^\m''^ l^^sm^sm((^-ni) , 

. ^ ( R r^-5R^\ . , , 

^'^[m^ 2R^ ) ~ ' 

A-) (R r=^ + 5i?V + 2i?3\ 

^0 2R^ ) cos^cos((/? - fit) , 



ci-) ( ^ , r^ + 5R'r + 2R^ \ 

= ^ [~M^ + ) sin^sin(99 - Qt) , 

and at r > it! it becomes 

d"*"^ -aflsin9sm((p — Q,t) , 

R 

= -2— //sin6'cos((^ - Qt) , 
^y = -2Riicos9cos{(p - nt) , 
^^+) = 2Rnsmesm{(p - fit) . 

VI. THE L = 1 METRIC PERTURBATIONS 

Once the actual gauge vector ^^"^^ is known, the MP are constructed via Eq. (31). Here we 
give the resultant expressions for /i^t\ first at the leading order (RH)^, and then at the next 
order {Rfiy (next subsection). Notice that all three time-space components hl^\h[f \hl^^ 
vanish at the leading order. The magnitude of these three components is oc Q. Notice also 
that hro and hr,p are smooth across r — R, and h^^ vanishes at both sides. 

A. Leading order (0°) metric perturbations 

Here are the values of h^^^ at the leading order {M/Rf = (RQ)^. 
Internal metric perturbations (2M « r < R) 

(-) f r R 

^ > = On i 

i?2 



^« ( ^ - j sin e cos{ip - Qt), 

hl;^ = 2// (^-^ + sin^cos((/? - Qt), 



hee^ = '^l^r^ - ^ ) sm9cos{(p - nt), 
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r 



R 

— — 2/x— sin^sin((^ — Q,t), 



hi'J = 0, 



'''tr — 'he ~ 'hip — ^■ 



External metric perturbations {R < r « 

R 

hlp = 4//— sin 6 cos((/? - Qt), 



— 2 fj,— cos 9 cos{(p — fit), 



R 

= -2/x— sin^sin((^ - fit), 



/,{+) _ /,{+) _ /,(+) _ /,(+) _ n 



'hr — 'he ~ 'hip — ^■ 



B. Next-order (Q^) metric perturbations 

Here we give the expressions of the three time-space components h^t^^ , h^^^ , h[^^ , which 
are all proportional to fl. All other MP components have vanishing contribution at this 
order, at both r < R and r > R. 

Internal metric perturbations (2M « r < R) 
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,(-) ^ ^ f 2r^ - WR^r + 5R^\ . , . , 
hfj. ' — z/iil I r„r>2 ) sin6'sin((^ — iU), 



hl^' — —2jAl ( ^-^^ ) cos^sin((^ — Vtt)^ 

'cos(</7 — Jli). 















(r^ + 107?V 







External metric perturbations {R <r 



/ii;r^ = -2/xl]i? y — J sin^sin((^ - f^t), 

— — J COS e sin(^ - nt) , 
/i^ ^ -2nnRi — — — J sin^cos(<^ - fit). 



C. Checking the above results 

The above expressions for h^^^ may be directly verified by applying the following checks: 

(a) The gauge condition h'^^ = is satisfied at both r > R and r < R; ^ 

(b) The vacuum field equation n/i^;, = is satisfied at both r > R and r < R; 

(c) The MP h^^^ is continuous at r = i?; 

(iv) The jump in ^^"^^^ agrees with Eqs. (32-34) above. 

Note that since we are only considering here the MP up to order Q,^, when implementing 
the checks (a,b) we may pretend that spacetime is fiat, i.e. M = 0. 

We should emphasize that the above checks (a-d) by themselves do not fully guarantee 
the validity of the derived MP, because they do not address the issue of correct asymptotic 
behavior at the EH. (In our construction, this issue is addressed by using the gauge vectors 
^j^^ with the correct asymptotic behavior at the EH and at infinity.) 



^The validity of this gauge condition at r = i? automatically follows from the continuity of the 
MP there, because the gauge condition only includes first-order derivatives. 
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VII. THE OVERALL METRIC PERTURBATION 



In this section we explore the structure of the overall harmonic-gauge MP, obtained by 
summing over all the tensor-harmonic modes. We shall consider here the strict weak-field 
limit, i.e. the leading order (M/i?)° = (Rfl)^ only. We then compare it to the trivial 
fiat-space MP solution. 

Consider first the fiat-space value of the even I — 1 harmonic-gauge MP. The particle is 
now assumed to be static (corresponding to Q = 0), located at 

r = R,9 = 7T/2,ip = ifiQ. 

A straightforward calculation of the flat-space / = 1 MP yields 

htt ^° = K'^^ = sin 6* cos(99 - (po), 

= 2//-^sin6lcos((/? - (po), 

3 

= 2//^(sin^)3cos((^ - ipo), 

for r < R, and 

hlP' = = 2/.^sin^cos(^ - (^o), 
hff^^^ — 2/iit! sin ^ cos(</7 — (po), 

h^+J'' = 2fiR{smefcos{ip-ipo), 

for r > R. The label "0" was introduced here to denote the flat-space quantities. All other 
components of vanish. 

Let us denote by A/i^j, the difference between the above fiat-space expressions for the 
even / = 1 mode, and the corresponding expressions of section VI for the Schwarzschild's 
weak-field limit. One finds (after substituting fit (fo) 
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Ahtt = -2n—sm9cos{ip - (po), (37) 
R 

Ahrr — 2/1— sin 9 cos{(f — (fio), (38) 

A.ho0 — —2ijlR sin 9 cos{(p — ^o), (39) 

A/i^<^ = -2iJ,R{sm9)^ cos{ip - ipo), (40) 
R 

Ahr0 — 2ijl— cos 9 cos((p — (fo), (41) 
r 

R 

Ahrip — —2//— sin 9 sin((^ — (^o) , (42) 
Ahe^ = A/it^ = Ahte = A/it^ = 0. (43) 



Note that although and Schwarzschild's hf^^y are non-smooth at r = i?, their difference 
Ah^y is smooth there. (This must be the case, because Ah^y is a homogeneous solution.) 

Consider now the sum over all modes. We shall assume here that apart from the even 
I — 1 mode considered here, all modes of the harmonic-gauge MP behave normally at 
the limit M — > 0. ^ (This assumption helps in providing a simple picture of the overall 
harmonic- gauge MP, though our main conclusions do not depend on it, as we discuss in 
the next section.) In the flat-space case, upon summation over all modes one recovers the 
standard fiat-space expression for the MP. We denote this overall flat-space harmonic-gauge 
MP by (here "0" stands for flat space, and "tot" for the total MP, i.e. the sum over all 

modes). In Cartesian coordinates x,y,z,t, for the flat background, takes the simple 

form 



^It appears likely (though it still needs be proved) that all radiative modes I > 2 will behave 
normally at the weak- field limit. For a circular orbit it can be shown that the / = mode is 
normal. It is likely that for a circular orbit the odd I = 1 mode is normal too. 
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with all mixed Cartesian components vanishing, where d is the distance between the particle 
and the evaluation point. 

Consider now the summation over modes in the Schwarzschild case (in the weak-field 
hmit). We denote this sum by Since the even I — 1 mode is assumed to be the only 

anomalous one, /i^*"*-* will deviate from just by the amount the Schwarzschild's even 

I — 1 mode deviates from its fiat-space counterpart, namely, 

/^r = /^?r^+AV- (45) 

We conclude that at the weak-field limit the overall harmonic-gauge MP is the sum of 
two terms: 

(i) The standard flat-space expression h^J^I;"^\ given in Eq. (44). This is a spherically- 
symmetric component, centered around the source particle, which decreases like 

(n) The anomalous component A/i^;^, given in Eqs. (37-43). This is an Z = 1 term cen- 
tered around the BH and oriented in the particle's direction. Note that /S.h^y is independent 
of the BH mass. 

Equation (45) is based on the above assumption that the even / = 1 mode is the only 
anomalous one. In case this is not the case, any additional anomalous mode will give its 
contribution to the anomalous component Ah^^. This is further discussed in the next section. 



VIII. DISCUSSION 

The term Ahf^,y in Eq. (45) demonstrates the anomalous behavior of the harmonic-gauge 
MP in the weak-field limit. This term is caused by the presence of the BH, but its magnitude 
is independent of the BH's mass M. Basically this term scales hke 

A V ~ ^ (46) 



24 



[see e.g. the tt or rr components, Eqs. (37,38)]. This demonstrates an interesting pathology: 
For fixed M and and fixed evaluation point in the neighborhood of the BH (though still 
in the weak-field region r >> M), the MP increases linearly with the distance R to the 
source (whereas naturally one would expect that the perturbation should decrease when its 
source is taken far away). 

Furthermore, for fixed M , n, and R, when one gets closer to the BH the anomalous term 
grows steadily as r~^. This growth continues as long as r » M. When approaching closer 
to the strong-field region, Ahf^u will approach a value of order 

AV ~ I^R/M^. (47) 

Consequently, even if the particle has a very small (but finite) mass /i << M, it will 
produce a large harmonic-gauge MP (A/i^^ >> 1) in the near weak-field neighborhood of 
the BH, if the particle is situated at a distance R » Rq, where Rq = M^//x. This may be 
achieved by either increasing R, or by decreasing M. 

In fact, if a finite-mass particle is situated at >> i?o it is no longer justified to consider 
the harmonic-gauge h^^y as a "linear perturbation", even if /j, << M, because Ah^i, >> 1. 
The entire perturbation approach thus breaks down in this case (provided that one insists 
on using the harmonic gauge!). Note, however, that the harmonic gauge can be used in a 
linear perturbation analysis for arbitrarily large R, if the mass n is taken to be infinitesimal 
(which corresponds to i?o ^ oo)- 

Detweiler and Poisson [3] recently found that the gravitational self force (in the harmonic 
gauge [8]) acting on a particle in a circular orbit around a Schwarzschild BH has a "Newto- 
nian" term, i.e. a term that does not vanish at the limit M —> (though its very existence 
does depend on the presence of the BH). Since the term in the overall harmonic-gauge 

MP is strictly isotropic around the particle (it is spherically symmetric), the "Newtonian 
self force" solely originates from the anomalous component Ah^,^. It is thus given by the 
contribution of Ah,j,i, to the effective "gravitational force" term 
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where SF denotes the perturbation in the connection, 




2^ ''(V;- + 



and we have used here the Newtonian hmit 5f . This only has r component in our case 
(all other components vanish as Q ^ 0). Substituting Eq. (37) in the standard expression 
for 5T one finds the self force 



This agrees with the analysis by Detweiler and Poisson [3]. 

The estimates in this section and the previous one were based on the assumption that 
the only anomalous mode is the even I — 1 mode. Even if this assumption is not satisfied, 
the pathological behavior described above still holds. Each additional anomalous mode will 
add a corresponding term to Ahf^u- Since each mode has a different angular dependence, 
it is not possible to cure the problematic terms (46) or (47) by adding other modes. For a 
generic angular direction, the small-r (but yet with r >> M) asymptotic behavior will either 
remain that of (46) — or perhaps become even worse [if additional anomalous modes exist 
which have a radial dependence more singular than (this appears very unlikely, however)] . 
Therefore, the pathological behavior of the harmonic gauge at the weak-field limit — e.g. the 
break-down of the linear analysis for R > Rq — is insensitive to the assumption made here. 

Next we comment on the domain of applicability of the expressions for the harmonic- 
gauge MP derived in sections VI and VII. Consider a circular orbit with some particular 
radius R » M (we assume here R « Rq, so the hnear analysis is apphcable). For too 
small r values, various terms proportional to M/r will become important. For too large r 
values, various retardation effects may become important. Therefore the domain of validity 
of the expressions derived in sections VI and VII isM<<r<<l/n. Note that in the limit 
M — > considered here, Q vanishes too, hence this range covers the entire r > domain. 

Throughout this paper we restricted attention to circular orbits. We expect, however, 
that the main qualitative results — the anomaly of the MP, the presence of " Newtonian self 



(48) 
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force" , and also the pathological behavior as R> Rq — will apply to other types of weak- field 
orbits as well. The explicit form of Ah^^, may depend on the orbit under consideration. In 
particular, for a non-circular orbit we do expect the anomalous component Ah^i, to include 
a contribution from Z = as well (however, the I — anomalous term is expected to be less 
singular than that of the even I — 1 mode, oc r~^) . The form of Ah^„ for non-circular orbits 
still need be explored. 
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APPENDIX A 

We exphcitly give here five of the six even, I — 1, exact static solutions for Eq. (5) 
in Schwarzschild spacetime. (The radial function of the last solution, will not be 

presented here as it is too long.) 

We first give the set of six modes with azimuthal number m — (i.e. the axisymmetric 
modes). The first pair of modes, denoted is given by 

with 

0(±) = ^(±)(r)cos^. 
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where 



A(-)(r)=r-M (A2) 

and 

^W(r) = (3/o) -2M + (M-OMl-2M/r) _ 
^ Next, ^^^^-^ are given by 

^(f^^ = 6lB^^\r)cose, (A3) 

where 

B^-\r) = r-2M (A4) 

and 

Finally, the pair takes the form 

^f±) = cW(r) cos^, 

ef±) = c(^)(r)sin^. 



The radial functions \ ■* are 



, -4M3 + 3Mr2 + 6r3 + 4MVln(r/M) 



6r 



^The factor 3/2, as well as the factor 3/4 in the expression for B^'^^ below, where introduced 
in order for these expressions to coincide with the corresponding expressions in section II in the 
weak-field limit. 
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and 




12M3 + 8MV - 3Mr2 + Qr^ + 8M^{M - r) ln(r/M) 
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The remaining two radial functions C^"*") , Cg are much more comphcate, and we shall not 
present them here. Both functions (which we have encoded in a MATHEMATICA file) are 
non-elementary, as they include the function Polylog[2, r/2M]. 

The other two sets of six modes corresponding to m = ±1 are obtained by straightforward 
rotations. One of these sets, the one which is mostly relevant to our analysis, is the following. 
The first pair i^^'^^-* is given by Eq. (Al), this time with 



{A^^^ as well as B'^^^ and Cjf"^ below are the same as specified above). The second pair is 



d^±) = 5*sW(r)sinesin(^. 



(A5) 



The last pair ^jf^^ now takes the form 



<:{C±) 



C^^^ (r) sin 9 cos (p , 




Cg (r) cos 9 cos ip , 



<:{C±) 



Cff ' (r) sin ^ sin (/? , 



,(C±) 



0. 



The last set of six m 7^ modes is obtained by transforming ip —>■ {ip + tt / 2) . 
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